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ABSTRACT 


This  wai  an  Investigation  of  the  mathematioal  methods  us^  in 
UDOFTT*  The  Ihveetlgatlon  has  pz<ovlded  an  evaluation,  both  frcn  a 
theoretical  and  an  experimental  aspedt,  of  tne  mathematical  methbds 
used  in  UDOFTT*  These  mathematical  Methods  Were  the  numerical 
methods  that  were  used  to  so]i.ve  the  simultaneous  system  of  nonlinear 
differential  equations  of  flight  for  the  F-IOOA*  The  investigation 
was  initiated  vlth  a  study  and  summary  of  the  mathematical  theory 
that  led  to  the  development  of  stability  charts  for  numerical  quadra¬ 
ture  fomulas*  Next,  the  theory,  uhloh  treated  a  single  llneax 
inhomogeneous  differential  equation,  was  extendsd  to  a  simultansous 
set  of  such  equations  and  then  to  a  primitively  nonlinear  differential 
equation.  This  last  equation  led  to  the  conjecture  that  stability 
charts  could  be  applied  to  nonlinear  differential  equations* in  order 
to  establish  both  stability  and  accuracy  of  their  solutions.  Finally, 
the  last  effort  experimentally  evaluated  the  accuracy  of  the  numerical 
methods  used  by  UDQFn,  as  applied  to  the  nonlinear  differential 
equations  of  motion  of  the  F-lOOA  trainer.  This  was  accomplished  by 
first  programming  certain  maneuvers  to  start  from  various  initial  con¬ 
ditions,  and  then  by  flying  the  trainer  throu^  the  progrmmned  maneuvers 
using  different  quadrature  intervals.  Quadrature  intervals  over  the 
range  of  100  msec,  to  10  msec,  were  used  and  the  resulting  flights 
were  compared  by  means  of  a  set  of  flight  varlablse  that  were  obtained 
In  digits,  print-outs* 
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FOREWCED 

The  mathematical  procedures  used  in  the  Universal  Digital  Operational 
Flight  Trainer  Tool  (UDOFTT)  have  been  under  investigation  since  the 
concept  of  real-time  digital  flight  simulation  was  first  envisioned.  The 
first  step  was  the  development  of  a  numerical  integration  procedure 
capable  of  solving,  in  real  time,  a  system  of  six  simultaneous  non-linear 
differential  equations  necessary  for  the  activation  of  a  flight  trainer. 

The  problem  was  investigated  through  the  theory  of  stability  charts  and 
their  use  in  uhe  synthesis  of  numerical  quadrature  formulae. 

Initially  the  theory  w^s  developed  to  treat  one  linear  inhomogeneous 
differential  equation.  The  natural  frequency  of  the  equation  was  computed 
and  for  this  frequency,  a  quadrature  interval  was  found  for  which  the 
solution  was  stable.  The  theory  was  then  applied  to  a  set  of  linear 
Inhomogeneous  equations  with  the  result  that  by  plotting  the  frequency 
curve  of  the  system  of  equations  on  the  stability  chart  an  immediate 
choice  of  quadrature  intervals  for  which  the  solution  will  be  stable  can 
be  made. 

The  present  study  was  an  experiments!  investigation  of  the  mathematical 
methods  used  in  UDOFTT  and  the  'variation  in  the  accuracy  and  stability  of 
these  methods  with  a  change  in  interval  length  and  under  the  stress  of 
different  flight  maneuvers.  A  by-product  of  this  investigation  was  the 
determination  that  stability  chart  theory  can  be  ^plied  to  a  primitively 
non-linear  differential  equation  when  there  is  close  agreement  between 
the  frequencies  of  the  true  solution  and  those  in  the  numerical  solution. 

The  flight  maneuvers  chosen  to  test  the  integration  procedures  "were 
of  both  high  and  low  frequency  types.  The  integration  interval  was 
varied  between  10  and  100  milliseconds  using  both  the  O^j  and  O^q  Cot 
quadrature  formulae.  The  results  showed  that  the  accuracy  and  stability 
for  all  intervals  chosen  using  bothfomulae  were  as  they  are  predicted 
by  the  stability  charts,  both  accurate  and  stable  for  the  10  to  70  milli¬ 
second  intervals. 

It  can  be  concluded  from  this  s'budy  that  the  O^o  Mod  Gurk  quadrature 
formula  can  be  utilized  for  the  simulation  of  most  airpljnes  even  though 
the  system  of  equations  describing  the  planes  are  not  "truly  linear.  The 
knowledge  of  the  accuracy  and  stability  of  the  chosen  quadrature  formula 
for  the  given  integration  interval  can  be  predicted  beforehand  through 
the  use  of  the  s'bability  charts.  Further  if  an  inci^ase  or  decrease  in 
computation  time  is  needed,  for  example  on  a  dual  or  triple  cockpit 
activation  system,  the  interv^  length  can  be  extended  to  70  milliseconds 
without  loss  of  accuracy  or  stability. 

— 

E.  M.  COLLETTI 

Con^uter  Branch 

U.  S,  Naval  Training  Device  Cen'ter 
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I,  INTRODUCTION 
1*  Historlcjll  Background 

Speculation  ea  to  the  feasibility  of  real-tine  ainralation  of 
airplane  flight  using  digital  conputers  dates  l^ack  to  the  middle  l^liO'e 
vhen  the  U«S.  Na'vy  Special  Devices  Center  laiinehed  a  atud^  to  examine 
the  problem*  At  first  glance  it  appeared  that  such  digital  Simulation 
should  be  feasible,  but  the  primitive  statue  of  the  digital  ooo^uter 
art  and  the  lack  of  certain  msthematical  assutancea  led  to  early  dis- 
eouragement,  at  least  for  the  immediate  future. 

In  19^0  the  fesaiblllty  study  was  re-instituted  at  the  Moore 
School  where  it  was  at  first  concluded  that  feasibility  was  marginal 
for  as  simple  an  airplane  as  the  F9F  but  that  feasibility  might  be 
accomplished  through  the  uiw'  of  certain  apeclal-purpoae  techniques 
aimed  at  increased  computational  speed*  These  techniques  proved 
successful  and  the  result  was  a  design  for  a  Universal  Digital 
Operational  Fll^t  Trainer  (UDOFT)  which  was  eventually  constructait  as 
a  trainer  tool  (UDOFTT ) * 

Concurrently,  the  mathematical  problem  underlying  successful 
digital  flight  simulation  was  solved  to  the  approximation  needed  at 
that  time*  Briefly,  the  mathematical  result  stated  that  the  stability 
Of  simulated  flight  could  be  guaranteed  to  be  essentially  the  same  as 
real  flight  for  all  possible  airplane  flight  conditions.  In  other 
words,  the  realism  of  the  flight  simulation  could  be  guaranteed  in 
advance,  regardless  of  the  flight  path  to  be  followed. 

With  the  deliveiy  of  UDOFTT  to  the  Navy  Training  Device  Center 
in  i960,  the  opportunity  was  presented  to  demonstrate  theee  results  in 
a  real  simulator .  The  present  study  has  undertaken  to  demonstrate  the 
validity  of  the  mathematical  proof  and  to  extend  the  result  to  an 
empirical  dsmonstration  of  the  accuracy  of  digital  fll^^t  simulation 
of  sn  FolDOA  alrplans  * 

2*  Technical  Background 

The  mathematical  proof  of  a  priori  stability  of  digital  flight 
simulation  reata  qpon  the  "quadrature  formula  stability  chart."  The 
latter  la  a  mapping  in  tha  conplex  plane  showing  tha  affect  of  digital 
computer  (or  "numerical")  solution  lulng  a  spaciflad  integrating  or 
"quadrature"  formula  on  a  eat  of  linear  ordinary  differential  aquatlona* 
More  specifically,  the  stability  chart  maps  tha  true  frequanclaa  of 
tha  analytic  solution  into  tha  correeponding  freqtianolas  of  the  nunarical 
solution.  The  stability  chart  diaplaya  both  tha  error  in  oBolllatlon 
period  and  the  error  in  daxplng  rate .  Formulas  exist  tdiich  permit  the 
oonputation  of  total  error  in  the  nunerlcal  eolation  of  Uneer  equations 
at  ary  spaciflad  time* 
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stability  la  aaeured  even  vhan  aocura<y  la  poor  aa  long  aa 
the  numerical  solution  ia  positively  (iajn)ed»  The  stability  ohart 
carries  a  curve  of  aero  damping,  whidi  is  therefore  the  boundary  of 
stability. 

Actually,  the  ft^uencies  referred  to  above  are  so-called 
"dljnenBionleoB  frequencies,"  i.e.  the  actual  frequencies  multiplied  by 
the  quadrature  Interval  or  time-step  by  uhich  the  .flight  solution  is 
advanced.  By  proper  selection  of  the  time -step  ai>d  quadrature  formula, 
it  is  always  possible  to  keep  the  numerical  solution  well  within  the 
stable  region  of  the  stability  ohart.  Unfortunately  thie  calls  for 
too  short  a  time-step  for  some  quadrature  formulae  and  existing 
digital  computers  could  not  keep  up  in  real  time.  Fortunately,  some 
quadrature  formulas  have  been  found  which  eatisiy  all  these  conditions 
but  at  some  snail  saorillce  in  accuracy. 

The  stability  chart  applies  only  tof  linear  ayeteme.  Airplane 
characteristics  are  only  qua eillnear,  i.e.,  they  change  but  the 
change  is  slow  enough  that  the  quadrature  formiila  treats  tliem  as  linear. 
The  results  of  this  study  indicate  that  the  aesumption  of  quasi" 
linearity  holds  for  the  P-IOOA  but  not  as  well  as  might  be  desired  by 
some  users  who  require  very  hi^  accuracy  as  well  as  stability.  For 
them,  the  solution  is  faster  digital  simulators  (or  mors  efficient 
programming  of  UDOFTT)  and/or  improved  quadrature  formulae.  Thus, 
although  the  present  study  has  confirmed  the  favorable  prediotions 
of  tbe  stability  chart,  further  effort  is  warranted  to  extend  the 
techniques  of  digital  fli^t  simulation. 
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n.  WORK  FERPOBMED 
la  Introduction 

Thlc  section  deeorlbea  the  work  perfozmed  for  the  evaluation  df  the 
mathamatlcal  methede  for  UD(^I  during  the  {wiet  year.  A  eunmary  of  previoue 
nathematlcal  studies  along  with  sons  recent  studies  are  presented  In  Section 
2«  Included  among  these  recent  studies  are  some  on  the  meaning  of  stability 
charts  for  nonlinear  differential  e<]uatlons.  Section  3  discusses  the 
experlaisntal  work  perfonaed  from  the  viewpoint  of  theory  and  method.  The 
maneuvers  studied  aie  listed,  the  Initial  conditions  md  the  control  positions 
for  the  maneuvers  are  discussed,  and  the  ooaputer  program  is  brought  to 
date* 

2.  Summary  of  Mathematical  Studjes 

The  first  Progress  Report  prepared  under  this  contract  was  devoted 
to  presenting  a  summary  of  the  previously  conducted  mathematical  studies  which 
led  to  the  development  of  stability  diarts  for  quadrature  formulas.  The 
suiaaary  deiaonstrated  the  applicability  of  the  stability  diarts  to  sets  of 
linear  hoDogenepus  differential  equations.  In  the  past  Hoore  School  reports, 
the  manner  In  which  forcing  functions  enter  into  t)w  numerical  solution  was 
described  only  for  eingLe  equations.  Tha  effect  of  each  forcing  functions 
''  on  sets  of  equations  was  derived  and  described  In  Progress  Rspoit  No.  2. 

Also  presented  In  Progress  Report  No.  2  was  a  section  on  Introductory  com- 
BSBtson  the  meaning  of  stability  charts  for  primitively  nonlinear  differen¬ 
tial  equatlona. 

The  following  seotions  present  a  brief  smamary  of  the  material 
referred  to  abovei 

2*1  Solution  of  Linear  Homogeneous  Dlfferentlel  Equations 

Thje  general  constant-parameter  linear  ayatam 

*l  •  •u’^  *  •••  *  Vn 


*n-‘nA*  —  *  V*n 

may  be  written  In  the  form 

I-  AX  (2) 

where  X  la  the  column  vector 
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•ad  A  la  the  aatrlx 


2«1*1  Solution  Using  Jordan  Cmonioal  Ibm 

• 

Equation  (2)  showa  that  the  deriratiTe  X  ia  cbtalnad  fron  the 
ordinate  vector  X  hy  a  linear  tranaformation  T,  which  In  the  X  CQordlnat* 
egrstea  la  repreeented  by  the  matrix  A.  When  the  proper  choice  of  ooordinat* 
•yatem  la  made,  the  niatrlx  repreaenting  T  takea  on  an  eapaolally  aliqd.a 
fomi  and  the  aolutlon  of  (2)  beoonee  quite  eaay. 

The  matrix  U  can  ba  obtained  ao  that 


X  -  UT  (3) 

I  -  U“^X  .  (U) 


Uaing  (3)  and  (U)  the  relationahip  between  A  and  B,  the  amtrloaa 
of  T  in  the  X  and  in  the  T  ayetema,  can  be  ooapated, 

I  «  IT^X 

i  -  U"^  •  IT^AX  •  u"\oi  -  BI,  (5) 

Then  the  following  theorem  waa  eatabllBhod: 

Lpt  X  and  T  be  related  b7  (3)  and  (U)*  Then  Z  la  a  ablution 
of  X  -  AX,  if  and  only  if  T  la  a  aolutlon  of  t  •  BI,  where 
A  and  B  are  related  by  ($) , 


The  objective,  thua,  la  to  find  a  coordinate  hyatea  in  vddch  the 
matrix  haa  a  minlaun  of  off-diagonal  eleewnta.  An  eapeoially  oonrenlent 
aituatlon  occura  when  T  haa  n  linearly  independent  eigenvectora 


with  elgenvaluea  1^, 

•••$  **  diagon^  elemanba, 

ipitM 


In  thla  oaae  matrix  B  ia  diagonal,  having 
and  It  ia  neoeaaary  only  to  ooaalder  the 


m 


Thaa*  the  genarai  aolutlon  of  (2)  In  thla  oaae  la  glTon  tgr 


(6) 
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(7) 


vfaar* 


VAien  the  charactoriatio  equation  has  multiple  roots,  the  Jordan 
oanonleal  form  reduces  the  solution  of  system  (2)  to  the  solution  of  systems 
of  the  foxm 


71  -Xyi 
^2  "  ^2  *  ^1 

4  •  %  Vl. 


Equations  (8)  are  of  such  a  form  that,  knowing  yj(t),  X^l(t) 
oan  he  obtained.  Thus, 


vhsre  o^  “  7^(0)  «)• 

For  additional  details  refer  to  Seotion  11*2  sad  to  i^ipendlx 
of  Progress  Report  No.  1,  M.S.R.  No.  60>23. 

2.1.2  The  Laple.oe  Trinsform  Point»-of-View 

Section  Il-b  and  Appendix  A»3  of  Progress  Report  No.  1,  H.S.R, 

No.  60-23,  show  that  Laplaoe>transfem  techniques  oan  also  be  used  to  sqIts 
equation  (2)  to  reach  the  sane  solution  as  that  obtained  in  Seotion  2.1.1. 

2.1.3  Nunerical  Solutions 

OiTsn  an  method  the  numerical  solution  of  X  •  AX  wss  found  to 

bs 


t  I 


5 


(10) 


lUplaelng  ty  OY^  and  multiplying  by  iT^  gave 
M  N 

Whan  A  la  diagonallaablaj  only  the  f ollovdng  ayBbem  needa  to  be 
ooaaldarodi 

»i  *  Vi 

Refer  to  Eq*  (6) 

4  "  Vei 


Thaaa  are  the  equationa  for  which  stability  charta  were  first  studied;  and 
80  it  waa  Icncwn  that 


y^(nh)  -  ^  Cjj  (e  }”  (1-1,  n)  (12) 

where  ^^h  are  the  dlmenaloiileas  frequenolea  whlob  the  nuurioal 

BWthod  aasooiates  with  the  dijunslOriLeas  frequencies 

For  the  eysten 
Yl  -  1  Yl 

y^2  •  ^  Y2  ♦  Yi 

.  Refer  to  Eq«  (8) 


ym  "  ^ 


the  aama  assumptions  ae  above  on  Xh  gave 
r  ^ik  ^ 

YjcC^-)  Pjk  (n)  (•  )“ 

whore  all  Pjj((n)  are  polynomlala  In  n  of  degree  mJ.  or  lower. 

Appendix  A-U  In  Progrooa  Report  No.  1,  M.a«R.  No.  60>23,  and 


(13) 
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Progress  Report  No.  2,  M.StR.  No.  6l*07|  dedorLbes  the  \iie  of  0^^  mod  Qurk 

to  solve  a  330  system  of  linear  homogeneous  differential  eq:aatloits>  aad.gl'veil 
the  0^^  mod  Qurk  quadrature  fomula  In  closed  f omii 

2.2  Numerical  Solution  of  linear  Inhomogeneous  Dlffeyentltl  Equatlona 

It  was  the  purpose  of  this  phase  of  the  work  to  study  the  nuiterl- 
cal  solution  of  linear  equations  with  time-dependent  forcing  functions  'tMch 
are  slowly  varying  with  a  view  to  determining  the  value  of  staBllity  charts 
In  predicting  the  stability  and  the  accuracy  of  such  solutions.  Prelimin¬ 
ary  results  Indicated  that  both  the  stability  and  the  relative  accuracy  of 
the  solution  of  inhomogeneous  equations  can  be  predicted  from  the  stability 
chart  I  and  whenever^  the  honogeneoua  solution  is  stable^  the  inhomogeneous 
will  also  be  stable. 

2.2*1  Solution  of  Single  Equations 

Application  of  a  given  open-  quadrature  formula  td  ths  hcno- 
geneouB  equation  a  ,  ^  ™ 

yielded 


B+p 


(15) 


idiere  s  »  max  (M,N),  ^  k  >  M  arui  the  b^  for 

k  >  N  are  set  equal  to  zero*  The  closed-fom  eolxitlon  to  (15)  is 

^  i  !»*»»  ^  ^ 


x(nh) 


(16) 


°k 


where  the  r^  are  the  roots  of  the  characteristic  polynomial  of  l.e.^ 
e 

r®  -  -  0,  (17) 

3"1 

and  the  are  constants  determined  from  the  s  initial  valuee. 


Application  of  one  of  the  open-closed  formulas  having  the  some 
stability  chart  as  the  0^  formula  to  ih  ■  lx  led  to  the  same  equation  (15) 

•when  the  proper  conditions  held  for  the  o.  (Refer  to  p.l-U)  of  M.S.R*  No. 
57-02).  “ 

*  Sin*  *e<p  •  *8-j+p  *  ^  '’a  (Hi) 
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To  determine  the  effect  of  4  forcing  function  on  the  solution^ 
via  mpplled  to 


X  ■  \x  ♦  f  (t)  • 

(18) 

This  yielded 

0 

*S4p 

■^1  “iVj.p 

*  ^9*p 

(19) 

where 

5  h  b.  f  f  (n-3)h)  1 . 

(Wa) 

In  cloeed-form  the  solution  is 

8 

"k  'k" 

*  k-1 

n  f  r 

•5  J 

(20) 

Note  that  if  Ir^^j  ■  2,  3»  ...|  »»  almoBt  the  entire 

tern  f^  is  applied  to  r^(orp^  while  the  other  frecpienoiee  ■*** 

are  only  slightly  excited.  Horeorer,  eren  these  small  excitations  das^ 


out  much  more  rapidly  than  the  r^  tern.  In  fact>  in  this  special  case 

It  follaiTB  that  the  nunerical  solution  of  the  inhonogeneous 
equation  (18)  has  the  following  dharacterlstiost 


(21) 


1.  The  solution  is  the  sun  of  two  ooinqponents. 

2.  One  ooa^nent  of  the  solution  (the  oomplenentary  term)  is 
identical  with  the  solution  of  the  oorrespoixting  homogeneous 
equation, 

3*  The  other  compohent  of  ths  solution  (the  particular  solution) 
is  composed  of  the  sum  of  the  sepiirate  contributions  Of  the 
forcing  function  at  eaoh  step  of  the  solution^  the  contribu¬ 
tion  at  each  step  being  ind^andent  of  all  the  other  eontrihu- 
tions  and  of  the  ooiqi^Iattentary  term. 

U«  At  the  n'th  step  of  the  solutionf  the  neff  contribution  of 

the  forcing  function  f^  depends  not  oxily  on  the  stability  chart 


bib  alao  on  tho  particular  quadratur*  forauLa.  usad.  This  ia 
the  only  difference  between  0|g|  and  Op^  fortalaa  with 

the  earn  stability  chart. 

Baoh  contribation  piropagattfs  exactly  as  the  ednpleaentary 

soltxtion.  All  of  the  frequenoiee  are  excited 

and  no  others.  In  particular ^  if  the  coa^pLesisntary  tana  li 
stable  in  the  sense  that  its  contribution  deoa3rs  asymptoti*. 
colly  to  sero  regardless  of  starting  conditions  x^f  3^,  ...i 

x^^t  then  the  particular  solution  is  also  stable  in  the 

sense  that  it  does  not  grow  without  bound.  Horeover, 
aoouraoy  late  in  the  solution  is  not  affected  by  errors  intro¬ 
duced  early  in  the  solution. 

6.  It  one  root  r^  greatly  exceeds  all  the  other  roots  in  absolute 

Mgnitudef  the  corresponding  frequency  will  be  the  only 

one  excited  to  any  appreciable  extent  and  alaost  the  full 
effect  of  the  foroli«  function  will  be  implied  to  this  exci¬ 
tation. 

7*  A  quadrature  foraula  with  a  good  stability  chart  will  hare 

X  and  jrjj  j  I’jjj  all  k  >  1.  The  quadrature  solu¬ 
tion  of  the  Inhonogeneous  equation  will,  therefore,  be  deter- 
nined  by  the  "natural  resonant  frequency**  of  the  corresponding 
bixaogeneous  equation. 

Ibr  additional  details  refer  to  Section  Il-$  of  Progress  Report 
Mo.  1,  H.S.R.  No.  60-23. 


2*2.2  Solution  of  a  SyEtwu  of  Eciuatione 


The  general  oonstant-coeffioient  inhonogeneous  aystea 


il  ■  *  •••  *  ‘lA  ‘  'i<») 

• 

4 

e 

5  ■  *  •"  * 

■agr  alsb  be  written  in  vectorial  foz»  as 


X  -  AX  +  F(t) 

where  X  is  the  oolunn  vector 


(22) 


(23) 
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F(t)  It  tht  column  vtotor 


■nd  A  It  tht  matrix 


to  equation  (23)  with  a  •  max  (M^N)  ^lelde 

8  t  a 

‘L  “iFn-k  "  "  r  ^  "n-k  ‘  ‘Zr 

k-1  W  k-1 

Any  m  dimensional  veotor  V  can  be  viritten  in  the 

(1) 

vhere  the  are  the  elgenveotoxe  of  an  arbltraxy  n  by  it  conatant  natrioe 
and  the  y^  are  linear  conbinatl'Ons  of  the  components  of  V«  Let  the  eigen- 

veotors  of  A  be  «•*!  ^re  is  atsooiated  vith  eigen¬ 

value 

Equation  (2li)  can  thus  be  urltten  for  the  firet  quadrature  ttep 


n 

X<1> 


D 


(  \ 

■'  ) 


0  XQ 

x(ih)  >2^  %(&  C(8-k)h] 
feo.  Vi 

*JI  \  ^ (H yj  ^ B-k)hix< J 

Wl  3-1 

(25) 

'V  ^  C(t-b)h3  . 
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AaBune  that  tlie  matrix  A  has  a»  dlatimbt  elganraLoea*  Then  there 
le  a  matrix  %  euch  that 


where 


ir^  AX-  -  D 

/\  0 

'  \ 


V 


Consider  the  tr.anaformation  X  •  X.  ^ia  known  to  be  the  matrix 
...  X^"^).  Hence  , 

fo\ 

.  1«>.  <>  .  •‘c. 


where  X 


Hence  equation  25  can  be  transformed  to 

8  ,  m 


,  x«)  . 

i\\ 

f 

9 

/ 

\ll 

,  etc. 


Keh)  -  (ah)  -  ^ 

^  \  J  I  (-W-J 


k" 


k-1 


(26) 


11 


«har« 


wter* 


[(B-k)hJ  [(i^k)hj  . 

Equation  (26)  may  ba  aaparatad  into  tha  oonponanta 

I(eh)  •  ^yj^(Bh),  ygCah),  y^Cihjj  , 

Iteration  than  prodicaa  that 
8 

^n-k)h^  , 


h  b,. 


k-1 


khara 


^j^n-k)hj  ■  (^j_|j»-k)hJ  , 

•  I^(n-k)hj  f  and  j  •  1,  »,,,  m 
Equation  (27)  can  be.  written  in  the  alnqiler  form 


vbere 


•jk  ”  “k  *  •>  ’'A 


and  -r- 


^  3n  ■  ^  ■ 


k-1 

The  aolutlon  of  equation  (26)  In  cloaed-form  !• 


y^C*!^)  •  ^  r 
^  k-; 


r  "  ♦  i  £  - 

>1^  k,! 


Hanoa  the  final  solution  to  aq;uatlon  (23)  la 
I(nh)  -  )Cl(nh)  --^yj^(nh),  ,  y,(nh))  ^ 


(27) 


(28) 


(29) 


(30) 
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Kltar*  tlM  7^(Ah)  «r«  dateinliMd  by  (29}« 

For  Additional  dotalle  ref«r  to  Motion  of  Progroso  ftoport 
Mo.  2.,  M.S.S.  Mo.  6147. 

For  tho  reauLti  of  applyl&gr  this  Mtbod  to  a  3  by  3  syaton  of 
oqiiationo  refer  to  Appendix  B  of  H.S.R.  Mo.  61-07. 

2.3  Introductory  COBunente  on  the  HeanliiK  of  Stability  Charts 
i[‘or  Mnd-tiyaiy  hoidlnear  ^fXerentlal  Eqantiona 

A  study  was  Inltlatod  to  detennlne  the  extent  to  which  the 
nathaaatloal  work  done  In  the  past  with  linear  differential  equations  sppllea 
to  nonlinear  differesbial  equations .  Work  was  perfomed  ttsl4ag  a  primitively 
Oonllnear  equation  and  a  slsqjle  nunerloal  awthod. 

The  nonlinear  equation  considered  was 

X  •  Ax  ♦  «x  .  (31) 


It  is  easily  werifled  that  the  analytic  solution  of  this  equation  Is 


(32) 


For  sufficiently  snail  e,  the  solution  can  be  approxiaated  as 

x(t)  ^  x^  j  •  ^^3) 

th  9ll  foznula  was  used  for  the  nunerloal  solution. 

The  n""  waOLtie  ofHt  as  given  by  ‘^ll  waa 

■  Vi  *  ^  ■w  “>  • 

‘  “nJ  *  “  •  <*> 

where  r  ^  l-^hX  ~e^  and  E  ^  he.  Note  that  the  definitions  of  r  and  pi 
are  idantically  those  that  appear  In  the  derivation  of  the  stability  charts* 

The  next  step  was  to  obtain  an  expreesion  for  the  solution  it 

any  point  In  tense  of  x_.  Based  on  equation  (jb)  and  the  assunptlon  that 

2  ® 

E  « 1,  it  was  proved  by  Induction  that 
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(35) 


For  E  <  4 
be  written  u 


p  X,  .  E  (yy^. . g 

1  and  with  e**^  eubstltuted  for  t,  equation  (35)  oaa 


x(nh) 


»«<» 


)inh 


1  ♦ 


he 

7^ 


(36) 


For  euffiolently  email  (ih  equation  (36)  le  approximated  by 

x(nl.)  ~  x„  fl  -  ./  (x  .  i^)]  .  Or) 

Comparlaon  of  equation  (37)  with  equation  (33)  flhowe  that,  with  the  reetrio- 
tiona  stated  above,  the  solution  of  a  flrs)i)>order  ellghtly  zonlinear  differ¬ 
ential  equation  obtained  by  the  numerical  method  0^  bears  an  interesting 
similarity  to  the  analytical  solution  of  the  eguatxon«  The  two  solutions 
differ  only  in  that  (33)  is  in  terms  of  X  and  (37)  is  in  terms  of  p. 


It  was  conjectured  that  stability  qharts  can  be  used  to  indicate 
the  conditions  for  accuracy  of  the  numerical  solution  by  noting  the  regions 
where  there  is  close  agreement  between  ph  and  Xh,  The  results  of  two 
oomputatlons  using  equations  (33)  and  (37)  are  given  to  illustrate  the 
accuracy  of  0^^  as  iqpplied  to  a  nonlinear  differential  equation. 

Example  1. 

x^  ■  1,  X  -  -1,  e  •  0,01,  p  •  0,1,  t  ■  1,  h  ■  0,1 

From  equation  (33)  x  ■  0.370 
From  equation  (37)  x  •  0,352 

Exampla  2. 

Xjj  ■  1,  X  -  -1,  e  «  0,01,  p  •  2,  t  •  1,  h  ■  0,1 

From  e(;piation  (33)  x  ■  0i370 
Prom  e(;piation  (37)  x  «  0.350 
By  0j2  (10  iterations)  x  •  0 , 350 

This  work  was,  of  course,  concerned  with  a  simple  equation  and  a 
Biiiq:ile  numerical  method,  ikiwever.  the  indication  is  that  more  oonqlex  altua- 
tions  should  be'  iixirestigated  in  the  hope  that  similarly  gratifying  results 
may  be  obtained. 


For  addl-Uonal  details  refer  to  Section  II-2  of  Progress  Report 
No,  2,,  M.S.R.  No.  61-07, 


lit 


8«k  CcBimentg  on  a  Runi-c-Kutta  Foapwula 

MotlTSted  ly  the  flz)dihg8  relating  to  0^^  aod  Ourk  com* 

putationa  prCeented  in  Final  Demonatration  Report  lS-2f  Jetluaxyi 
19$9t  "An  Eacporijnentol  Analog-Digital  Flight  Simulator,*  a  preliminary 
study  was  conducted  relative  to  the  significance  the  stability- 
chart  criteria  as  applied  to  the  secure <y  of  solutions  of  nonlinear 
diffarentlsl  equations «  The  following  Runce^Kutta  type  formula  was 
first  considered t 

(36) 


idiere 


^*1 

*1^  ‘^o 


ii  .  ,  ill  .  , 
*1  *l 


il 


1  h  X  ^ 
7  "  *1 


_  iii  .  1  ^  •  ii 
*1 


Iv 


i  h  *1 


iii 


This  formula  was  uaed  to  obtain  the  solution  for  the 
linear  equation 


Xx. 


The  true  solution  of  this  equation  ia  known  to  be 


x(t)  “  e 


Xt 


^0  ® 


.Xhn 


*0  ® 


.an 


09) 

(Uo) 


irtiere  z  4  Xh  ,  h  "At,  and  t  ■  hn. 

The  solution  using  formula  (38)  is  obtained  as  foUowsi 


■  >“0 


hXx. 


ii 


iii 


1  ^  •  i  .  1  V,,,  i  ,  1  ,2  , 
7  ^1  7  hXxj^  ■  ^  B  X 


1 


iv 


1 V..;  ill 
1 


«|hx. 


Jhloc, 


3T 


0 


ill  .  1  X 
1  "lu**o* 
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*1  •  *.  *  “o  *  ii  *0  *  ^  *c  *  ir  *c 


0  0 

Lit  \  ••  *1.  Thtn  (39)  beooM0 
i  •  -tt, 

and  with  •  1  ai«l  h  ■  1  (lil)  glaat 

U 

-  0.37$0. 

Tha  tma  solution  given  tgr  (bO)  la 
X  (1)  -  a*^ 

•  0.3679* 

The  percent  error  la  1.9J(. 

Next  the  ncnllnaar  aquation^ 

i  •  -X  ♦  i  x^, 

was  oonaldarad. 


(bl) 


a2) 


The  tma  solution  for  (b2)  with  x.  •  1  la 

o 

X  (t)  -  ~  . 

l-fs^ 

Trcm  the  Taylor  series  aapanalon  for  f(x«t)f  it  foUowi 
that  for  the  first  interval 


BenoS)  for  equation  (b2)  with  x^  *  1  aid  h  ■  1 


(b3) 

(bb) 


X  ■  -1  +  ^  Initially 

X'^*0*$  at  the  and  of  tha  Intarial 


16 


Uelng  the  Runge-Kutta  type  formule  (38)  to  obtain  the 
solution  for  (I42) 

-  hi^  •  -1  +  I  (1)^  -  -0.5000 
•  I  hx^^  -  I  [.5000  ♦  I  (-.5ogo)^l  -  00125 
-  I  hxj^^  -  ^  £-3125  ♦  j  (.3125^  -  -0.0879 
*1^’'  *  *  f  I  (-.0879)^  -  OJ0229 


Xj^  -  1  -.5000  +  .3125  -JO079  ♦  />229 
-  0.7U75 

The  true  solution  for  t«l  as  gliien  bT  (U3)  ie 

X  (1)  -  0.5379. 

The  percent  error  is  39$. 

Since  Xh*0  initially,  a  more  aocurate  solution  night  have  been  expected 
baaed  on  the  Runge-Kutta  stability  chart.  Apparently  the  large  Change 
in  Xh  over  the  interval,  1  .e .,  the  strong  nonlinearity,  causes  serious 
error  in  this  Runge-Kutta  type  formula* 

The  classical  Runge-Kutta  formula  for  solving  an  equation 
of  the  form  x  ■  f(x)  ie 

“  *0  ♦  A.x  (W) 

where  ^x  -  ^  ♦  2k2  ♦  2kj  ♦  k^) 

with 

kj^  -  f  (x^)h 

kj  -  f  h 

k,  -  f  (x„  -  h 

k^  -  f  (*0  ♦  kj)  h. 


I 
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Far  equation  (U2)  thla  gava 

^  .5000 

kg  -  -OJt687$ 
kj  -  -OJ^Ta^A 
-  -0.388356 

^  X  •  -0  Jt6l8zi 

•  l-04i6l8 
-0.5382. 

The  percent  error  la  O.D69C. 

The  daaaloal  Runge-Kutta  gives  a  nudi  better  result  than  does  farmula 
(38)  «hen  applied  to  equation  (U2)  even  though  the  two  fomulas  have 
the  aame  stability  chart.  The  reason  for  the  difference  is  not  as 
yet  fully  understood. 

For  conparlson  the  mod  Ourk  quadrature  foraols  was 
used  to  obtain  a  solution  for  vU2)  over  the  saoie  Interval.  It  gave 

•  0  Jt667 

When  the  proper  past  values  were  used.  The  percent  error  Is  13^^ 
vdileh  Is  hl^er  than  what  might  have  been  ejqpected  based  on  the  0^^ 

SAd  Qurk  stability  chart. 

The  results  of  this  preliminary  study  point  out  the  need 
for  determining  the  significance  of  the  stability-ohart  criteria  idien 
applied  to  the  accuracy  of  the  solution  of  nonlinear  differential 
equations. 
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3*  BbcpcriJnental  WbA 
3»1  Theory  and  Method 

The  quadrature  methods  used  in  Mere  studied  and 

evaluated  hy  comparing  results  obtained  froiri  programmed  autbpllot 
memeuvers  using  vari.ous  quadratui^  intervals*  Several  flight  mdneu- 
vers  were  selected  for  studj*  and  for  demonstration  of  the  propertlee 
of  digital  flight  simulation.  On  the  one  liacd,  maneuvers  Mere  selected 
which  would  demonstrate  the  advantages  of  digital  elmlationi  and  on 
the  other  hand«  maneuvers  were  sought  vhieh  would  display  the  llaltations 
of  digital  simtilation* 

Each  maneuver  was  started  from  one  of  four  ln>alr  positions* 
Each  position  was  defined  by  a  set  of  ihitlal  conditions j  l.e.^  the 
infoimation  needed  to  initialize  flight.  Autopilot  programs  were 
generated  to  provide  repeatable  control  poaitioiin  as  a  function  of  tine* 
Each  maneuver  was  lloHn  with  several  different  quadrature  interval^ 
using  the  0^^  mod  Qurk  and  the  Gu'k  quadrature  foimiles*  The 

quadrature  intervals  used  were  10  msec.,  20  msec.,  $0  msec.,  70  msec*, 
and  100  msec.  (It  should  be  noted  that  the  duration  of  the  computation 
oyele  was  held  constant  while  the  quadrature  jntervals  were  changed* 
While  this  increased  the  computation  time  Vith  the  shorter  quadratura 
intervals,  it  had  no  effect  on  the  results.) 

In  order  to  use  the  above  quadrature  intervals,  certain  quan¬ 
tities  which  depend  upon  the  quadrature  interval  had  to  be  changed  for 
each  interval.  These  quantities  included  the  0^,  mpd  Oark  coefficients 

d,  e,  and  f,  and  the  O^qC^  mod  Qurk  coefficient  d^*  (A  list  of  the 

constants  with  their  values  for  the  various  quadrature  intervals  can  be 
found  in  Appendix  E  presented  in  Progress  Report  No.  3^  MSR  lo.  61-13 
and  this  report.) 

While  the  maneuver  was  being  run,  analog-output  recordings 
were  made  of  the  following  quantities i  u,  v,  w,  p,  q^,  1^,  n^,  and  H 

(an  exception  to  this  occunrod  during  the  short-period-osclllation-in- 
yaw  maneuver  where  r  replaced  H  and  replaced  w) «  During  the  flight 
these  quantities  were  monitored,  and  at  its  end  the  minimum  and  the 
maximum  values  of  the  quantities  wore  printed  out  to  provide  a  scale 
for  the  tracings.  In  addition,  values  for  the  quantities  to  be  studied 
were  periodically  stored  in  the  number  memory  and  then  printed  out  in  a 
block  at  the  eiid  of  the  maneuver.  The  fraquenoy  with  vhLoh  the  values 
were  stored  depended  upon  the  maneuver* 
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3*2  DiBflusslon  of  I'fork  and  Aoproiofa 
3*2*1  M«neuv»r8  Studied 


1 


The  maneuWs  that  «er«  atudltd  fall  into  thr«e  groups  ss 

follows! 

These  maneuvers  Involve  sesentislly  steady^stats 

a.  Steady  level  flight 

b.  Steady  rate  of  cUnib. 

Group  2t  These  maneuvers  involve  twoHliaenelonal  transient 
moiion . 

a.  Phugoid  oscillation. 

b*  Short  period  oscillations  in  yaw  introduced  by  applying 
a  pulse  to  the  rudder. 

c.  Short  period  oscillations  in  pitch  intx*oduoad  ly  applying 
a  pulse  to  the  elevators. 

d .  Snap  roll . 

e .  Stall 

Proto  These  maneuvers  involve  three>dijaenaional  motion, 
a.  Immelman  turn. 


b.  Split  S. 

The  maneuvers  were  started  from  "in-flight*  initial  conditional 
The  four  initial  conditions  requited  were  classified  by  the  altitudes 
and  the  Indicated  air  speeds  ss  follows! 

Initial  Condition  I  30^000  feet 

36$  knots 

Initial  Condition  II  30,000  feet 

200  knots 

Initial  Condition  III  20,000  feet 

1$0  knots 

Initial  Condition  IV  10,000  feet 

500  knots 

The  awneuvers  classified  by  initial  conditions  are 
I.  1.  Steac^y  level  flight. 

2.  Short  period  oseiUstions  in  yaw. 
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3  •  Short  period  oaclUstlons  in  pitoh 
k  *  Sn«p  HoU . 

5.  Phugoid  otolUation. 

II .  1 .  SpUt  S 

III.  1.  Stall. 

IV.  1.  Stea4^  rate  of  ollad). 

2.  Lmelman  turn. 

3.2.2  Obtaining  Initial  Conditions  for  the  ysneuvers 

In  order  to  evaluate  the  numerical  methods  uaed|  it  was 
necaeaary  to  start  each  maneuver  from  identical  conditions  for  each 
flight  vdth  the  different  quadrature  Intervals  so  that  they  could  be 
realistically  compared .  If  the  initial  conditions  had  not  been 
identical  for  each  execution  of  a  maneuveri  there  could  be  little 
useful  comparison  of  the  results. 

A  set  of  initial  conditions  is  the  information  needad  to 
oontlnua  computations  from  a  given  point  in  time .  Bach  of  the  f cur 
seta  of  initial  conditions  was  obtained  as  follows t 

A  pilot  took  off  in  the  trainer.  The  discrete -input  switches^ 
which  could  not  be  set  while  the  trainer  was  on  the  ground |  were  set  to 
their  pre-determined  positions  As  soon  as  all  the  discrete -input 
switches  were  prcperily  set  the  trainer  was  flown  to  >*iat  seemed,  from 
Infomstlon  observed  on  the  instruments,  to  be  a  steady  level  flight 
at  the  specified  altitude  and  Indicated  air  i^eed.  Vlien  this  condition 
was  attained,  the  controls  were  released.  Beosuse  there  are  no  con¬ 
trol  forces  in  the  P-IOOA  trainer,  releasing  the  controls  was  equivalent 
to  placing  constant  control  positions  in  the  memory  and  using  them 
to  fly  the  trainer.  The  flifht  was  continued  with  these  constant 
control  positions  for  about  one  minute .  During  this  test  fli^t, 
the  instruments  (xt  the  instructor's  console  were  observed  to 
determine  thether  a  steady  level  flight  hed  been  attained.  If  it  had 
not,  the  pilot  took  over  control  and  tried  to  bring  the  trainer  to  a 
better  position}  and  the  above  test  was  repeated.  If  steady  Isvel  flight 
had  been  attained,  the  FRBRZE  switch  was  turned  on  and  a  "print-out"  pro¬ 
gram  was  read  into  the  instruction  memory.  Tl^s  print-out  program  con¬ 
sisted  of  CLA  instructions  oalllng  for  the  contents  of  the  nuaber-isemoxT 


♦tiscrete-input  switches  were  set  so  that  the  flights  would  be  flown 
with  a  "clean"  airplane  with  low  weight. 

Refer  to  Appendix  C  in  Progress  Report  No.  3|  MSR.  No.  61-13. for  a 
list  of  the  proper  switch  positions.  ’ 


loc«tlons  that  are  of  interest.  The  cortputer  waa  then  put  In  CONT 
PRIMT,  and  the  quantities  were  printed  out.  These  quantities  were 
then  studied  more  closely,  and  the  quantities  considered  necessary 
to  initialiae  a  flight  frwa  an  in-air  position  were  punched  on  in¬ 
put  caixis  for  UDOFTT.  The  chosen  quantities  were  then  tested  by 
using  them  to  Initialise  a  steady  level  fli^it  with  the  constant  con¬ 
trol  positions  that  were  obtained  tttn  the  print-out . 

3 .2 .3  Obtaining  Control  Positions  for  Maneuvers 

The  control  positions  for  the  autopilot  program  for  the 
various  maneuvers  were  obtained  in  the  following  ways 

A  pilot  flow  the  trainer  to  a  steady  level  flight  corresponding 
to  the  proper  initial  condition.  I'Jhen  the  steady-level-f light  con¬ 
dition  had  been  attained,  FREEZE  was  turned  on  and  the  proper  set  of 
initial  conditions  was  reed  into  the  number  memoxy.  The  controls  were 
then  set  to  the  proper  positions  to  agree  with  the  values  in  the  initial 
conditions.  The  analog-output  recorder  was  started,  FREnZE  was  tumad 
off,  and  the  maneuver  was  executed.  JPhile  the  maneuver  was  being 
executed,  the  control-position  inputs  were  monitored  by  the  program  and 
then  multiplexed  out  to  the  analog-output  recorder .  The  analog-out- 
put  recorder  produced  tracings  of  the  control  positions,  gnd  the 
min  .-max .  monitor  program  provided  the  minimum  and  maximum  values 
of  the  recorded  functions  from  which  scales  were  determined.  The 
control  positions  were  then  approximated  by  straight  linos  and  thus 
mode  available  for  use  in  the  autopilot  program.  (A  88ii^)le  approximation 
is  shot-m  in  Appendix  D  in  Progress  Report  Ko.  3i  MSR  No.  61-13; 

3.2 Test  Program 

3 .2  fU  Introduction 


The  test  program  was.vritten  to  fly  the  trainer  trm  an 
“in  air"  condition  through  a  predetermined  maneuver  at  aty  one  of 
several  quadrature  intervals  using  one  of  two  quadrature  formulas 
while  monitoring  several  of  the  variables  in  the  flight  equations. 

In  order  to  do  this,  the  autopilot  program  provides  the  control 
positions  as  a  function  of  time  and  provides  analog  and  digitel 
outputs  for  the  monitored  variables.  It  would  have  been  desirable 
to  have  one  complete  program  to  perform  these  tasks  fcr  each  maneuver 
at  each  quadrature  internal.  However,  ccnsideraticn  of  the  large 
number  of  cards  vftjioh  would  be  needed  for  these  programs  prompted 
the  division  of  the  programs  into  four  groups  which  could  be  used 
in  many  combinations  to  form  the  large  number  of  programs  needed* 

The  four  parte  of  the  autopilot  program  were  celled t 

1 .  Initial  Conditions 

2.  Basic  Program 

3*  Maneuver  Programs 

U<  Corrections  for  Quadrature  Intervals 
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The  autopilot  program  used  irt  the  inVe atigption  Is  listed, 
in  Appendix  G  which  starts  in  Progress  Report  No.  3>M5R  Ko4  61-13# 
and  is  continued  in  Interim  Technical  Reoort  MSR  61-19 »  Ihis  appendix 
contains  the  complete  listing  of  the  Basit  Program,  the  Maneuver 
Program  (except  for  the  scaling  values)  tor  each  of  the  nine  maneuvers, 
end  two  examples  of  the  Corrections 'for-Quadraturc -Interval  programs. 
These  programs  are  listed  in  preliminary  and  in  final  code .  The 
preliminary  code  'is  £iven  in  terms  of  the  mnemonic  code  for  the 
instructions  end  the  symbol  for  the  nianber -memory  address.  A  discussion 
of  the  instructions  and  their  mnemonic  codes  can  be  found  in  PROORAIMINO 
Hi^NUAL  FOR  TIS  UIXDFTT  C0MPDTP.R.  Aupist  19$9.  Appendix  F,  which  was 
started  in  Progress  Ifteport  l^o.  3#  nSR  61-13#  end  is  ocntinued  in 
Interim  Technical  Report,  MSR  61-19,  contains  a  list  of  the  symbols 
which  have  been  introduced  in  these  programs  and  their  definitions 
are  given  in  terms  of  the  variable  represented  by  the  symbol#  however 
the  symbol  is  also  used  to  designate  the  number-memory  location  Wiere 
this  variable  is  stored.  Other  symbols  used  are  takan  from  UDQjFTT 
SIMULATION  PROGr^M,  FINAL  REPORT.  MAY  I960  VOLUMES  I  and  II. 

The  autopilot  program  is  entered  from  the  permute  portion 
of  the  F-IOOA  program,  and  the  F-IOOA  program  is  re-entered  via  the 
usual  permute  exit  which  has  been  placed  at  the  end  .of  the  autopilot 
program.  It  was  necessary  to  eliminate  aome  of  the  F-IOOA  program  to 
provide  instruction-memory  locations  for  the  autopilot  program. 

Because  the  maneuvers  are  flown  in  the  air  with  the  No-Puel-Ueplation 
Look  ON,  the  Land-  Air-  Crash  DedeionB  and  the  Mesa  of  Fuel  portions 
of  the  F-IOOA  programs  were  removed  dth  no  rceulting  effect  on  the 
coIl^)utatlons  * 

3 .2 .1| .2  Initial  Conditions 


Group  1  of  the  autopilot  program,  Initial  Conditions,  contains 
four  sets  of  flight  variables  which  are  used  to  initialize  the  aim- 
ulator  to  the  flight  condition  from  which  the  maneuver  is  to  start. 

Bach  set  contains  those  variables  needed  to  determine  a  given  in¬ 
flight  condition.  A  discussion  of  these  conditions  and  the  methods 
for  obtaining  them  can  be  found  in  Section  3.2  J.  and  3.2  42. 

3 .2  .Ij  .3  Basic  Program 

Group  2  of  the  autopilot  program,  the  Basle  Program#  contains 
those  parts  of  the  autopilot  program  which  are  common  to  all  of  the 
maneuvers  to  be  performed.  Thus  group  2  contains  one  program  which  1« 
used  for  all  maneuvers.  This  program  provides  the  analog  outputs,  the 
digital  outputs,  and  the  timer,  as  well  as  various  routing  and  control 
instructions . 

The  autopilot  program  provides  eight  analog  ou^uts  for  the 
analog  recorder.  During  eaeh  computation  cycle  eight  variables  were 
scaled  and  multiplexed  out.  The  scaling  waa  performed  to  obtain  a 
magnification  of  the  variation  of  the  variables.  To  obtain  the  proper 
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■eAlng  a  trial  maiksuver  waa  executadj  and  the  digital  outputs  frcn  thli 
■anevnrer  provided  the  naxliimm  and  the  minlnun  values  that  the  variables 
attained  during  the  aanouver.  The  malbwum  and  the  minlniim  values  for 
.each  variable  were  then  averaged  to  flpd  the  nidpoint  of  the  variation 
for  that  variable*  the  difference  between  the  majdmum  and  the  nlnlnum 
provided  the  Magnitude  of  the  variation  whloh  was  used  to  determine 
the  anount  of  shifting  needed  to  bring  this  variation  to  the  most 
significant  portion  of  the  output.  The  midpoints  add  drifts  were  then 
added  to  the  Maneuver  Program  so  that  in  all  subsequent  flights  the 
midpoint  for  each  variable  waa  subtracted  from  the  variable  and  the 
difference  was  then  shifted  left  to  a  more  signlfieent  position  and 
noltiplexed  out. 

The  digital  outputs,  which  were  provided  by  the  autopilot 
program,  weret  1)  a  Hat  of  maximum  values  attained  by  each  of  eight 
variables,  2)  the  times  associated  with  these  nuDdnums,  3)  the  min¬ 
imum  values  attained,  U)  the  times  associated  with  the  minimums,  and 
5)  a  list  of  values  of  these  variables.  The  maxlmuns,  the  miniaums, 
their  associated  times  were  provided  by  the  Maximum-Mininum  re¬ 
corder  in  the  Baisic  Progrm.  During  each  computation  cycle  the  Max- 
Ifflun-Miniinum  recorder  compared  the  present  value  of  each  variable  with 
the  number  stored  as  the  maximum  value  for  that  variable,  and  the 
larger  of  the  two  was  stored  as  the  new  maximum.  If  the  present 
value  was  as  large  as,  or  larger  than,  the  past  maxinpom,  the  present 
time  was  stored  as  the  time  that  the  maximum  occurred,  replacing  the 
time  previously  recorded.  Similar  comparisona  were  made  to  detendne 
the  mdnlinum  values  and  their  associated  times. 

The  list  of  values  was  obtained  by  periodically  storing  the 
values  of  the  variables  in  a  block  of  memory.  An  index  and  the  Tally 
Register  were  used  to  inorease  the  storage  addresses  in  forming  the 
list.  A  prinb-out  program  waa  provided  to  enable  the  removal  of  this 
list  and  the  lista  from  the  Maximum-Miniimai  recorder. 

A  timer  was  needed  in  the  autopilot  program  heoauae  the 
control  positions  were  obtained  as  functions  of  time.  The  timer  used 
10  msec*  (scaled  BlU)  as  the  unit  of  time  (i.e.,  1  in  the  lUth  binary 
position  from  the  left  was  eqpal  to  10  msec.)  for  all  maneuvers  except 
the  Phugoid,  where  20  msec,  was  used.  During  each  oamputation  cycle 
the  timer  added  the  quadrature  interval  being  used  to  the  time  for  the 
previous  cycle  and  stored  the  result  as  the  time  for  the  present  eyola. 
The  end  of  the  desired  flight  waa  determined  by  oom{>aring  the  time 
oaloulsted  in  this  maraier  with  a  predetexmined  number.  When  the  flight 
was  over  the  autopilot  program  sat  the  crash  flag  so  that  the  trainer 
entered  the  mode  of  oamputation. 

The  Basic  Program  also  containsd  olear-add  orders  which 
replaced  the  mnltiplex-in  orders  for  the  control  positions  in  the  F-IOOA 
pngrmrn  The  olesr-add  orders  called  for  the  control  positions  from 
Mmory  locations  rather  than  from  the  cockpit  oonbrolsj  this  enabled 
the  autopilot  program  to  provide  the  oohtxol  positions  nsoessary  to 
perform  a  maneuver* 
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Maneuver  Pro^ramg 


Group  3  of  the  autopilot  program,  the  Maneuver  Programs, 
contains  a  program  for  each  of  the  nine  maneuvers  to  be  performed# 

For  a  list  of  these  maneuvers  refer  to  Section  3«2«1«  Each  of  the 
programs  contains  a  list  of  the  breakpoints,  the  slopes,'  and  the 
intercepts  for  the  control  functions  used  to  give  that  maneuver.  The 
method  used  to  obtain  these  functione  i^  described  in  Section 
The  control  positions  are  calculated  in  the  Function  Generav>r  with 
time  as  the  independent  variable*  They  are  then  scaled  khid  stored  for 
later  use  by  the  F-IOOA  program.  As  a  result  of  the  scaling  usSd  the 
control  positions  were  stored  with  the  same  precision,  10  bits,  if 
they  were  provided  via  the  analog«-to-digital  convolvers  used  In  the 
trainer# 


The  Maneuver  Programs  also  contain  the  list  of  midpoints  and 
shift-left  orders  for  the  analog-output  scaling.  Although  these  orders 
are  contained  in  the  Basic  Program,  their  magnitudes  were  set  according 
to  the  maneuver  being  performed*  They  were,  therefore,  stored  with 
each  maneuver* 


3#2*1|*5  Corrections  for  Quadrature  Intervals 

Group  U  of  the  autopilot  program,  Corrections  for  Quadrature 
Intervals,  contains  a  program  for  each  maneuver  a't  each  quadrature 
interval  used,  or  fifty-four  programs  in  all.  Each  program  is  a  small 
group  of  nurabera  and  instructions  which  had  to  be  changed  when  the 
q^iadrature  interval  was  choiged.  These  include  constants  which 
involve  the  quadrature  interval  such  as  tha  weighting  coefficients  in 
the  integration  formula,  the  time  increment  for  the  timer,  and  shift 
orders  for  scaling  made  necessary  by  changes  in  these  constants*  For 
a  list  of  the  quadrature  intervals  used  and  a  discussion  of  the  con¬ 
stants  which  depend  on  them  refer  to  Section 


3.2.1|*6  Oy^  Flight  Program 


The  test  autopilot  program' described  above  was  based  on  and  used 
with  the  F-IOOA  simulation  program  which  used  the  0.^  mod  Gurk  quadrature 
fomula*  The  additional  program  used  during  the  study  was  a  modification 
for  the  F-IOOA  simulation  program.  This  modification,  when  used  in  con¬ 
junction  with  the  F-1O0a  program,  provided  a  simulation  program  using  the 


jicd  Qurk  quadrature  formula. 


The  classical  integration  formula  is  formed  by  the  oyclio 

combination  of  the  classlWl'^u^prediction  fomula  and  the  classical  C— 
closure  formula.  In  thie  mannw  each  new  value  of  the  variable  is  estinated 
by  the  0 ^formula,  the  corresponding  derivative  is  computed,  and  the 
improved'^ stimate  of  the  variable  is  found  using  the  formula.  The 
0„C,,  mod  Gurk  quadrature  fomula,  which  was  used  for'^e  study,  is  an 
l^rovenient  over  the  classical  0^.0.^.,  The  mod  Qurk  formula  i* 

given  in  detail  in  Apendix  E.  ^  ^ 
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To  obtain  an  F-IOOA  simulation  program  that  used  the 

mod  Ourk  quadrature  formulaj  modifications  were  made  on  the  existing 
simulation  program  which  uses  the  Ourk  quadrature  formula. 

The  0^^  integration  routine  was  replaced  by  the  routines.  The 

program  was  then  passed  through  twlOSi  the  flret  time  ueing  0^  and  the 

second  time  using  for  each  cycle.  SinOe  two  sets  of 

values  for  each  variable  were  involved  in  each  cycle  ,  i.e.j  the  first 
and  second  estimateSi  new  variable  storage  locations  and  provisions  for 
determining  vMoh  vklue  was  needed  were  added.  It  was  found  that  all 
the  ooaputationa  needed  for  the  mansuvers  could  be  run  in  this  way  In 
the  50- nsec.  Interval  provided. 
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in.  RESULTS 


In  this  section  the  results  obtained  for  the  quantities 
of  primary  Interest  for  each  msneu^r  td.ll  be  discussed.  A  Hat 
of  these  quantities  and  their  symbols  afe  glvtn  below: 


u: 

velocity  vector  along  the  airplane  X-exls 

w: 

velocity  vector  along  tha  airplane  Z-axls 

ps 

roll  rate 

pitch  rate 

cosine  of  the  angle  between  the  Xrsxls  of  the 
plane  and  the  fixed  Z-axls  of  the  esrth 

air- 

cosine  of  the  angle  between  the  Z-axls  of  the 
and  the  fixed  Z-axls  of  the  earth 

airplsns 

ft 

Sideslip  angle 

H: 

cltitude  above  grotmd 

The  maneuvers  were  flown  using  the  0^^  mod  Gurk  and  °30Sl 

mod  Ourk  quadrature  formulas  td.th  quadrature  Intervals  of  10,  20,  ^0, 

70,  and  lOO  msec.  The  results  thus  obtained  are  presented  In  drawings  (in 
Appendix  lO  tdtereeadi  drawing  shovis  the  results  with  one  of  the  two 
quadrature  formulas  and  the  siQiroprlatc  quadrature  Intervals .  When 
there  Is  little  difference  between  the  results  obtained  with  the 
various  quadrature  Intervals,  some  of  the  results  are  aidtted  from 
the  drawing  in  order  to  preserve  clarity. 

The  comparison  of  the  results  for  a  given  maneuver  with  one 
of  the  qusdratxire  formulas  Is  based  upon  the  magnitude  of  the  difference 
between  the  results  obtained  with  each  qiiadrature  interval.  The  per¬ 
cent  error,  or  percent  difference,  figures  are  not  given  since  they 
depend  to  a  large  degree  upon  the  magnitude  of  the  solution  aa  wall 
ae  upon  the  magnitude  of  the  error.  The  frequency  of  the  solution  in 
those  maneuvers  which  involve  oscillations,  such  as  the  osolUatlons 
In  pitch,  is  noted  so  that  comparisons  can  be  made  between  the  various 
frequencies  obtained  and  the  frequencies  predicted  fan^  the  stability - 
cliert  criteria . 

Tha  comparisons  between  the  results  with  the  0^^  quadrature 

formula  and  the  quadrature  formula  are  also  based  on  the  magnitude 

of  the  differences  In  the  results.  For  these  conperlsona  the  results 
for  each  of  several  selected  maneuvers  at  the  ^O-msec .  quadrature  Inter¬ 
val  are  given  on  the  same  Figure  for  direct  comparison  of  the  alm- 
ulatlona.  For  some  maneuvers  the  difference  curve  for  these  two  Em¬ 
ulations  is  also  shown.  The  maneuvers  for  which  such  figures  are  given 
were  selected  to  indicate  the  type  and  magnitude  of  the  differences 
encountered  . 
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1*  Stetdy  Flight 

For  the  Ste«<ly-Flight  m^heuvei'  u  Mi#  aelacted  is  the  vsrlsble 
of  priinsry  Interest;  therefore  evaluation  df  the  data  obtained  was 
directed  toward  u.  The  results  are  displayed  in  Fibres  1*2 » 

The  steady  flight  maneuver  was  obtained  fay  simulating  flight 
for  '$0  seconds  with  constant  control  positions.  Digital  sanpllngs 
were  taken  at  one  second  intervals  during  the  flight.  The  0^^  formula 

was  used  with  two  quadrature  intervale,  50  and  100  msec.  The  data 
frcm  the  lOO-mseo.  run  Indicates  inaccurate  results;  however  the  solution 
was  marginally  stable  and  the  inaccuracies  are  not  sufficient  to  rule 
it  out  as  a  possible  simulation.  Thedita  from  the  50-msec,  run  In¬ 
dicates  a  stable,  acceptable  flight.  There  were  indications  that 
underflow  was  affecting  the  results  at  50  msec.;  thus  the  smaller 

intervals  were  not  used.  The  0,_C-,  formula  was  used  with  the  five 
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intervals,  and  there  was  veiy  good  agreement  among  the  results  obtained 
with  all  intervals  •  The  meximura  deviation  between  the  results  with 
ID  mstc.  and  100  msec,  is  lese  than  1  in  the  fourth  octal  place,  and 
the  results  v.  ■*h  the  other  intervals  fall  between  these  two.  The 
effects  of  unc  i._'^low  on  the  results  is  quite  apparent,  particularly 
at  h  -  10  msec,  where  no  variation  is  present.  The  difference  between 
the  results  obtained  with  the  tvro  quadrature  formulas  is  very  small y 
except  for  those  at  100  msec,  where  the  0.qC,.  eimulstion  gave  con¬ 
siderably  better  results.  ^ 

2.  Short  Period  Osclllationa  in  Yaw 


For  the  Oecilletlona-  In-Yaw  maneuver  was  selected  as  tha 
variable  of  primary  interest;  therefore  evaluation  of  the  data  obtained 
was  directed  toward  V'  .  It  should  be  noted  that  4*  itself  is  not  ob¬ 
tained  directly  by  0^^  (or  O^qC^^)  computations;  but  it  is  obtained 

from  the  ratio  kv/u,  vdiere  k  is  a  constant  and  u  and  v  arc  obtained 
directly  fay  O..  (0»C,^)  comoutatlona .  The  resulte  are  displayed 
in  Figures  3-7.  ^ 

The  maneuver  was  performed  fay  aoplying  a  triangular  pulse  to 
the  centered  rudder-control  position  at  time  t*5  seconds  during  a  steady 
level  fll^t  at  30,000  feet  and  Mach  0.98.  The  flight  duration  waa 
50  seconds,  but  digital  samplings  were  obtained  only  over  the  time 
Interval  t-l^.2  to  lli  seconds.  Two  types  of  pulses,  one  of  moderate 
and  the  other  of  large  amplitude,  were  used.  Karly  studies  indioated 
that  the  maximum  amplitude  of  f  obtained  with  the  moderate  pulsa  waa 
considerably  less  then  that  reported  in  the  MET  report.  Final  Dim- 
onstration  Raport  Ij5-2,  January,  1959,  "An  Experimental  Analog- 
Dif^ital  Flight  Simulator.".  Since  approximating  the  MIT  condition  was 
of  interest,  the  Oscillstions-in-Yaw  maneuver  using  the  large  pulse 
was  added  to  the  test  maneuvers. 


The  moderate  pulse  had  a  duration  of  1  sedond  and  ah  amplitude 
of  approximately  l/h  of  the  allovjable  dtttlectlon  fixm  the  cantered  rud¬ 
der-control  position.  It  produced  a  damped  oscillation  in  4>wtth  an 
average  period,  of  approximately  2.5  seconds  and  a  damping  factor  of 
about  0,2,  giving  a  n  of  approximately  -0*2+2. 5ii  'fhi  sollitions  with 
both  quadrature  formulas  at  all  quadrature  intervals  Were  stable  and 
similar.  The  principal  Variation  among  the  results  is  a  slight  ini- 
crease  in  the  period  with  decreasing  amplitude}  this  is  most  apparent 
vdth  the  10-msec. -interval  data.  It  was  suspected  that  underflow  and 
round  off  were  responsible  for  the  change  in  period  for  the  simulations 
with  h  ■  10  msec.,  so  difference  methods  were  used  to  approximate  the  ef** 
foots  of  round  off  for  the  0^^  simulation.  Because  jji  ■  kv/u  and  u  is 

essentially  constant  during  this  flight,  the  estimation  was  based  on  thd 
integration  of  V,  The  portion  of  the  curve  studied,  i.e.,  in  the  vicinity 
of  t  -  11,8  sec.,  was  selected  because  it  is  a  part  of  the  curve  with  a 
relatively  large  derivative  and  would  provide  a  con¬ 
servative  estimate  for  the  round-off  error.  It  was  found  that  for  this 
portion  of  the  flight  approjilmately  $%  of  the  contribution  of  the  contrl** 
Ixibion  of  the  derivative  to  the  integral  was  lost  because  of  round  off* 

A  much  higher  percentage  vras  lost  near  the  top  of  the  curves  where  the 
deidvatives  were  smaller.  The  presence  of  underflow  is  clearly  shown 
by  the  flattening  of  the  final  peak  of  the  0^^  slrmlation  with  h  ■  10 

msec.  The  absence  of  change  in  period  for  the  maneirver  with  the  lar¬ 
ger  pulse  also  tends  to  confirm  that  the  problem  is  related  to  underflow. 

The  large  pulse  had  a  duration  of  2  seconds  and  an  amplitude 
which  was  approximately  equal  to  the  allowable  deflection  from  the  cen- 
tei’ed  rudder-control  position.  It  produced  a  damped  oscillation  in 
with  an  average  period  of  approximately  2*ii  seconds  and  a  damping  factor 
of  about  0,26,  giving  a  p,  of  about  -0»26+2.6i,  The  solutions  with  both 
quadrature  fotmuilas  at  all  quadi'ature  intervals  were  stable  and  similar. 
Thore  is  such  close  agreement  among  the  results  with  the  different 
quadrature  intervals  that  only  the  results  with  10,  50,  and  100  msec, 
are  shown.  This  close  agreement  indicates  that  the  simulation  with  the 
50-insec.  quadrature  interval  is  fairly  accurate  and  certainly  meets  the 
requirements  for  a  satisfactory  simulation. 

For  a  given  quadrature  interval  the  0^^  simulattojr  and-  the 

O30  ^31  close  agreement  for  the  oscillation  in  yaw 

excited  by  the  moderate  pulse  and  by  the  large  pulse.  The  two  simula¬ 
tions  for  and  their  difference  vith  -toe  large  pulse  and  the  50-msec. 
quadrature  interval  are  ahawa  in  Figure  ?•  The  difference  curve  in(}l- 
cates  that  the  two  formulas  gtors  a  slightly  different  phase  response  -to 
the  forcing  function  bxit  that  the  simulations  of  the  unforced  portions 
of  the  oscillation  differed  only  in  the  magnitude  of  the  oscillations* 

The  scales  for  Figure  ?  show  that  the  amplitude  difference  for  the  t-wo 
simulations  was  about  $%  with  0^^  having  the  larger  magnitude. 
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Ab  a  whole  the  resvilte  obtained  are  in  close  agreertent*  No« 
where  was  difficulty  of  the  nature  reported  by  HIT  encountered,  and  no 
explanation  can  be  deduced  for  their  difficulty#  It  is  noted  that  for 
our  maneuver  the  moderate  rudder  pulse  produced  a  inax  Ax  de¬ 

grees  and  the  large  pulse  produced  a  max/^7  ;^i2  degrOeb,  while  that  for 
the  group  solution  in  the  MTT  report  is  al^ut‘''17  degrees.  It  is  also 
noted  that  17  degrees  exceeds  the  limit  of  15#9  degrees  placed  j  Uj/ 
in  the  P-IOOA  program#  '  #  / 


) 
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3.  Short  Period  OeclllatloPB  In  Pitch 


For  the  OBoillatlons-ln-Pltch  oaaevtver  q-  vu  Selected  ai  t:h« 
variable  of  primary  Intereet;  therefore  evaluatibn  of  the  datta  obtained 
vae  directed  tovmrd  .  The  reeults  aire  displayed  In  Plguree  8-11 

The  maneuver  vas  performed  by  applying  a  triangular  phlie  at 
time  t  ■  5  seconds  to  that  elevator  position  required  to  maintain  steady- 
level  flight  at  30,000  feet  and  Mach  O.98.  The  pillse  had  a  dtiratlon  of  1 
second  and  an  anqplltude  of  (>-. 0200000 )g*  .  The  flight  4tlril.'tion  2k 
seconds,  but  digital  samplings  vere  obtained  only  over  the  tlmd  interval 
t  X  4.2  to  9.1  seconds.  The  pulse  produced  an  osclllatlcm  In  q.  with  an 
average  period  of  approxliaately  1  second,  and  a  4  of  approxima-Mly 

nD.7  +  6.5i 

The  results  with  the  0--  formula  are  displayed  in  Figures  8-9. 

The  results  with  100  msec,  show  a^lltudes  which  are  much  greater  than  those 
with  the  other  quadratture  intervals;  however  the  envelope  of  the  curve  in¬ 
dicates  that  the  oscillation  is  growing  in  a  linear  manner  instead  of  the 
exponential  growth  encountered  with  an  unstable  solution.  This,  together 
with  the  fact  that  the  minimum  and  maxlniuB  values  occurred  at  approximately 
11-12  seconds,  indicates  that  although  the  solution  is  unsatisfactory  because 
of  the  large  oscillations,  it  is  stable.  The  results  with  10,  20,  ^0,  and  70 
msec,  are  displayed  in  Figure  0  and  are  similar  to  each  other, 

The  results  with  the  O^o^Sl  using  10,  50>  100  msec,  are 

displayed  in  Figure  10.  The  three  tracings  are  similar;  however,  as  with 

there  Is  a  slight  decrease  in  the  period  as  the  quadrature  interval  is  in¬ 
creased.  This  variation  is  expected  in  part  from  stability-chart  criteria, 
but  it  seems  to  be  present  to  a  greater  degree  here  because  of  the  nonlinearity 
of  the  F-lOQA  flight  equations.  There  is  also  more  daaqplng  with  the 

than  there  is  with  the  O^^.  However,  the  most  notable  variation  from  the  0^^ 

results  is  that  the  lOO-msec.  results  track  the  others  quite  well  with  the 
OjqC^]^  while  with  the  0^^  the  lOO-msec.  results  are  wild.  This  difference 

seems  to  indicate  that  tte  O^qC^^^  formula  gives  better  solutions  to  the 

slightly  nonlinear  equations  when  the  larger  values  of  uh  are  reached. 

If  \  for  the  oscillations  in  q^  is  detexmlned  by  assuming  qlose 
agreement  between  the  true  solution  and  the  numerloal  solution,  l.e., 
between  \h  and  uh,  with  h  «  50  msec.,  a  \  of  approximately  -O.7  +  6.51  is 
obtained.  The  locus  of  \h  for  all  h  on  the  O^g-mod-Gurk  stability  chart 


^Subscript  6  denotes  octal  representation  of  numbers. 
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.!(  -0.7  +  6.51), 


100  nusieo.  Figure  1  (in  this  Section)  shows  the  O^j-nodMIurk 


is  a  straight  line  from  the  origin  throng  the  point  \h 
l.e.,  for  h 

stability  chart  (and  the  0  C  -mod-Gurk  stability  chart  since  It  Is  the 

same  sb  that  for  the  wltn  this  line  dravn  and  \h  tor  h  «  100  nteec., 

70  neec.,  and  50  msec,  marked.  It  1b  c^uite  apparent  that  the  fOr  h  ■ 
100  msec*  is  close  enoiu^  to  the  branch  contour  for  inaccurate  reaults  to 
be  expected  =  For  h  *  70  msec.,  however,  the  Xh  is  6tifficlentiy  far  a^y 
from  the  branch  contour  for  relatively  good  results  to  be  expected.  Ifdte 
clIso  that  an  increase  in  freq^uency  with  Increase  Ih  h  is  ihdlcated  by  ihs 
chart.  Such  a  shift  is  present  in  the  res\ilts  obtained;  hOweVet*  for  the 
0.^  sjlmulatlon  it  is  slightly  greater  than  that  indicated  by  the  stability 
cmrt.  The  stability  chart  predicts  a  solution  with  a  of  approximately 
-0.12  +  .751,  SDd  Figure  9  indicates  that  the  imaginary  part  of  the  \ih 
for  the  actiial  siraulatlon  using  0^.^  at  h  =»  100  msec*  was  approximately 
0.83.  However,  the  agreement  betwen  the  predicted  and  the  actual 
frequency  for  the  70-msec,  simulation  is  much  closer,  i.e.,  0.47  ve*  0.474. 
This  difference  may  be  caused  by  the  effects  of  nonlinearity  becoming 
slgnificeuit  near  the  branch  contour.  The  solution  with  the  0^0^31 

wild  for  h  »  100  msec.,  and  the  frequency  obtained  from  the  graph  is  qjrf.te 
close  to  that  predicted  by  the  stability  chart.  This  indicates  that  the 
O30C31  simulation  is  not  effected  as  much  by  the  nonllnearltles  of  the 

system  as  is  the  0^^  simulation. 


The  results  with  50  msec*  for  both  of  the  quadrature  fotmulas 
are  displayed  in  Figure  11.  The  main  variations  are  that  there  is  more 
damping  with  the  there  is  a  slight  increase  In  the  period  with 

the  toward  the  end. 
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4.  Snap  Roll 

For  the  Snap'-Roll  maneuver  p  and  n^  vere  selected  as  the 

variablea  of  primary  interest;  therefore  evaluation  of  th6  data  obtained 
vas  directed  toward  p  and  .  The  results  are  displayed  in  Figures  12-15. 

It  took  about  5  seconds  to  perform  the  roll,  after  which  the 
8imu3.ator  was  allowed  to  continue  flight  with  the  initial  Control  positions 
for  about  15  seconds.  During  the  **6teady  fli^t"  portion,  the  results  for 
both  p  and  n^  with  the  various  quadrature  intervals  ahd  both  quadrature 

formulas  are  so  nearly  equal  that  only  the  50-®fl®c.  points  are  shown .  Dur¬ 
ing  the  roll  portion  of  the  meneuver,  the  results  are  also  ih  close  egfeenent 
—  with  tracking  a  little  better  for  lu  with  the  particular*- 

ly  so  with  the  100-rasec.  interval  where  it  is  again  noted  thet  the 
seems  to  give  slightly  better  results  through  a  violent  maneuver.  ^  ^ 

5.  Fhugold  Oscillation 

The  phugoid  oscillation  Is  one  In  which  there  Is  a  large  amplitude 
vsurlation  of  air  speed,  pitch,  and  altitude,  but  d\irlng  which  the  variation 
of  angle  of  attack  is  slight.  Because  the  periods  are  relatively  long,  the 
effects  of  inertia  forces  and  damping  forces  eire  small.  The  phugoid  can  be 
thought  of  as  a  slow  Interchange  of  kinetic  and  potential  energy  about  some 
equilibrium  energy  level,  i.e.,  the  phugoid  can  he  considered  &e  a  long- 
period  oscillation  of  the  airplane's  velocity  as  it  attempts  to  re-establish 
the  equilibrium  condition  from  which  it  has  been  disturbed.  The  period  of 
the  phugoid  is  sufficiently  long  and  the  response  of  the  airplane  is  such 
that  negative  damping  of  this  Oscillation  has  little  bearing  on  the  pilot's 
opinion  of  the  flying  qualities  of  the  airplane.  For  this  reason  the 
military  often  has  no  requirements  on  the  damping  of  the  phugoid,  and  as 
a  result  many  military  airplamo  will  exhibit  unsteble  phugoid  characteristics 
within  some  range  of  their  operating  conditions. 

Included  with  the  various  maneuvers  run  during  earlier  phases  of 
the  project  was  one  for  phugoid  oscillation.  The  phugoid  maneuvers  vere 
performed  by  incrementing  the  integral  of  the  rate  of  climb  or  the  air  speed 
a  few  seconds  after  the  start  of  a  steady  fli^^t  using  Initial  Condition  I 
and  allowing  the  flight  to  continue  for  a  few  minutes.  With  small  changes 
the  trainer  quickly  resumed  a  steady  flight  (but  not  necessarily  at  the 
original  altitude).  This  effect  was  apparently  the  result  of  underflow  in 
the  calcvilation  of  u  and  H.  When  the  change  was  large  enough  to  produce 
an  oscillation,  the  trainer  went  Into  an  unstable  phugoid  oscillation. 

Figure  l6  indicates  the  types  of  results  obtained.  A  stable  phugoid  was  net 


*  The  discussion  in  this  paragraph  is  based  on  material  in 
Perkins,  C.D.,  and  Hags,  R.E.,  Airplane  .Performance 
Stability  and  Control.  New  York,  Nw  York:  John  Wiley 
wd  ^ons,  Inc .  /  1949  (50),  Chapter  10. 
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obtalmd  with  this  Initial  condition;  however  as  Was  pointed  out  above 
the  presence  of  an  unstable  pbugold  for  some  flight  conditions  of  a 
Bllltary  airplane  is  not  unusual,  so  it  was  assumed  that  these  tesults 
were  cooslatent  with  the  characteristics  of  the  airplane. 

There  are  sevexel  changes  in  the  airplane's  pareaeters  that 
can  be  made  to  Increase  the  dancing  of  a  phugoid.  These  Include  increas¬ 
ing  the  mass,  increasing  the  parasitic  drag,  or  moving  the  center  of 
gravity  to  a  more  favorable  position.  Therefore,  as  an  endeavor  to  obtain 
a  stable  phugoid  oscillation  the  phugoid  tests  were  run  from  a  new  condi¬ 
tion  which  Incorporated  the  above  changes  in  the  alrplaxie  pemameters.  To 
insure  high  drag  the  condition  was  taken  at  a  moderate  altitude,  12,000 
feet,  with  the  tip  tanks  on  and  the  speed  brakes  out.  The  addition  of  the 
tip  tanks,  which  were  full,  and  an  increase  in  the  fuel  load  provided 
an  Increue  in  the  mass  of  the  airplane  and  a  shift  in  the  center  of 
gravity.  The  condition  was  obtained  at  full  throttle  with  the  afterburner 
Off  which  produced  a  flight  at  Mach  0.5. 

For  this  test  the  phugoid  was  Initiated  by  decreasing  u  by  itO 
knots  at  time  t  <■  l4  seconds  during  a  steady  flight  using  the  phugoid 
condition  described  above.  The  dxiratlon  of  the  flight  was  322  seconds 
during  which  several  of  the  flight  variables  were  sampled  at  seven-second 
intervale.  Of  primary  interest  were  u  and  H,  \dilch  were  selected  as 
the  variables  that  would  best  exhibit  the  characteristics  of  the  phugoid 
oscillation.  The  results,  as  shown  in  Figures  17A  and  17B,  indicate 
the  ^presence  of  a  stable  phugoid  oscillation  superimposed  on  a  gradual 
climb.  The  period  of  the  oscillation  was  approximately  one  minute,  and 
the  damping  was  positive  and  appreciable. 

The  results  for  H  indicate  a  large  amount  of  disagreement 
between  the  simulations  with  vsurlous  quadrature  Intervals.  If  the  re¬ 
sults  for  u  are  considered,  however,  we  see  that  there  are  increasing 
amounts  of  round  off  leering  to  underflow,  indicated  by  the  flattening 
of  the  cxirves,  as  the  quadrature  intervals  were  reduced.  It  is  also  of 
interest  that  there  is  close  time-agreement  between  the  points  irtmre  the 
H  curves  diverge  and  the  points  where  the  u  curves  beccme  dominated  by 
vuaderflow.  It  was  decided  to  run  further  teste  to  substantiate  the 
statement  that  the  poor  agreement  in  the  simulations  was  caused  by  the 
presence  of  underflow  rather  than  the  characteristics  of  the  quadrature 
formula  used. 

In  an  attenpt  to  reduce  the  effects  of  underflow  modifications 

were  made  for  the  cODputation  of  the  derivatives,  u,  w,  1^,  m^,  n^,  and 

the  rate  of  climb  and  the  integration  of  these  quantities.  These  modifica¬ 
tions  involved  the  change  of  scale  factors  and  the  use  of  two  words  per 
variable  in  order  to  increase  precision.  The  Increase  in  precisian  for 


these  conqputatldtiB  ranged  from  tvo  td  six  hits  vlth  a  slx>hlt  Increase 
for  u  and  H. 

The  phvigold  tests  described  above  were  then  run  vlth  these 
changes  In  the  program.  The  results  obtained^  Figures  17C  and  17D  show 
very  good  agreement  for  all  quadrature  Intervsas;  and  there  Is  also 
close  agreement  between  these  results  and  the  results  vlth  h  «  XOO  msec, 
without  the  Increased  precision.  This  shows  quite  clearly  that  the 
difficulties  encountered  in  obtaining  phugoida  and  that  the  foor 
agreement  among  the  simulations  at  various  quadrature  Intervals  were 
caused  by  the  effects  of  round  off  and  underflOWi  Tbe  good  fesultS 
obtained  with  Increased  precision  and  the  phtigoid  fkaeuyef  demonsttates 
the  ability  of  the  0  mod  Qurk  formula  to  handle  lOW-ftequbncy  oSeiUao 
tlons  of  the  type  encountered  in  the  Phugold-OsciUation  maheuver. 

The  0„C_.  mod  Qurk  formula  was  not  used  for  the  Fhugoid- 
Oscillation  maneuver  since  it  had  already  been  shown,  for  exan^e  in 
the  Steady-Bate -of  >Cllnib  laaneuveri  that  for  nonviolent  maneuvers  the 
mod  Qurk  and  the  0^^  Qurk  formulas  give  very  similar  re¬ 
sults.  Thus,  the  inclusion  of  an  O^o^^Sl  maneuver  would 

not  have  provided  any  information  ttot  had  not  already  been  obtained. 

The  main  objective  was  to  investigate  the  behavior  of  the  0^^  mod  Qurk 

formula  during  the  simulation  of  a  maneuver  vtaoee  dimensionless 
frequency  falls  near  the  origin  and  slightly  to  the  left  of  the  x  ■  0 
line  on  the  stability  chart. 


tor  the  Solit-S  maneuver  w  and  p  were  selected  as  the  varlablea 
6t  prtikary  Interest)  therefore  evaluation  of  the  data  obtsihed  wea 
directed  toward  W  and  p*  The  results  are  displaced  in  Figures  16*22. 

For  w  with  0^^  the  results  M.th  10,  20,  and  $0  mate,  are  in 

close  atireemsnt,  and  the  lOO-msec.  results  tratk  these  falrljr  well. 

With  the  the  agreement  is  sasewhat  better.  With  the  50«mseb. 

quadrature  interval  the  results  with  the  two  formulas  are  nsarljr 
equal,  as  can  be  seen  in  Figure  20.  The  larger  variations  at  the  two 
places,  6*7  seconds  end  18*20  seconds,  are  probably  caused  by  the 
discrete  saenllng  of  ringing  functions  with  sli^t  phase  shifts. 

For  p  all  of  the  results  with  both  formulas  are  in  very 
close  sfreement. 

7.  StaU 


For  the  Stall  maneuver  u  and  L,  were  selected  as  the  variables 
of  primary  interest)  therefore  evaluation  of  the  data  obtained  was 
directed  toward  u  and  1^.  The  results  are  displayed  in  Figures  23*27. 

The  elevator  function  which  was  applied  to  the  simulator  to 
produce  the  stall  was  analogous  to  continuously  pulling  back  on  ths 
"stick"  in  ths  trsiner  i*!lle  flying  steady  fllg-ht  at  20,000  feet  and 
Mach  0.33*  Stall  occurred  at  approximately  8  aeconda  after  starting  the 
maneuver,  but  the  flight  was  oontlnuad  until  t-^  seconds. 

The  results  for  u  with  the  Oy^  formula  and  quadrature  inter* 

vals  of  10,  20,  $0,  70,  and  100  msec,  are  in  close  agreement  for  the 
first  18  seconds.  After  that  the  results  begin  to  diverge.  The 

results  are  quite  simllsr  to  those  fron  the  0^^  except  thet  the  pattern 

of  divergence  is  different*  Moat  notable  is  the  dccraaaa  in  the  value 
with  the  lO-maec.  interval  toward  the  end*  As  can  ba  seen  from  Figure  2$ 
the  ^0*m6ec.  results  for  both  formulss  sz*e  in  quite  close  sgreement. 

The  results  for  1^  with  the  0^^  formula  and  all  five  quad* 

rature  intervals  are  in  close  agreement  during  the  first  18  aeconda, 
and  then  they  befln  to  diverge.  The  0^  Si  results  are  quite  similar 

to  'Uie  Ojj  ones  exceot  that  the  pattern  of  divergence  is  different* 

With  the  ^OHBsec.  quadrature  interval  the  results  with  both  formulss 
are  in  quite  close  agreement. 


Since  the  point  of  dirargence  of  the  reeults  seechs  to  coireB** 
pond  to  the  point  of  minlmpi  velocity,  it  is  felt  that  the  variations 
are  caused  by  the  effects  of  anall  differences  in  the  flights  as  the 
plane  slides  into  a  dive  with  the  nose  up. 

8,  Steady  Rate  of  Climb 


For  tl:e  Steady«Rate»of-C3iiab  maneuver  1^  and  H  were  selected  as 

the  variables  of  primary  interest)  therefore  evaluation  of  the  data  ob¬ 
tained  was  directed  toward  1^  and  H.  The  results  are  displayed  in 

Figures  28-31. 


The  results  for  1^  with  quadrature  intervals  of  10,  20,  and 

$0  msec',  are  displayed  in  the  figures.  The  TOnaSec.  results  trhok  the 
JO-mseo.  ones  so  closely  that  they  are  not  shown.  The  lOO-msec.  results 
oscillate  and  go  positive,  but  the  minimum  amd  maximum  values  indicate 
that  the  results  are  stable.  While  there  is  close  agreement  among  the 
displayed  results,  the  results  with  a  smaller  interval  slightly  lead  those 
with  a  larger  one.  There  is  dose  agreement  between  the  results  with 
the  two  formulas. 


The  results  for  H  with  quadrature  intervals  of  10,  20,  ^0,  and 
100-4n8ec,  are  displayed  in  the  figures)  the  TOnnsec.  reeults  track  the 
50-nsec.  ones.  There  is  close  agreement  among  the  various  results  ex¬ 
cept  for  those  with  lOO-msec.  With  the  lOO-msec.  quadrature  interval 
the  simulator  went  into  a  dive  instead  of  climbing  and  crashed  at  about 
20  seconds  with  the  0^^  and  at  about  38  seconds  with  the  This 

difference  in  crash-times  is  the  principal  variation  between  the  results 
with  the  two  quadrature  fomulas.  This  again  shows  tendencies  for  better 
behavior  with  the 

The  start  of  the  erratic  behavior  of  1^  for  the  simulation  with 

h  ■  100  msec,  was  coincident  >tth  the  application  of  the  stick  motion  to 
the  simulation.  As  has  been  mentioned  before,  the  oscillations  in  q^^ 

and  1^  encountered  during  the  tests  have  frequencies  such  that  the  simu¬ 
lation  with  h  ■  100  msec,  sends  \h  close  to  the  branch  contour  on  the 
032-mod-Qurk  stability  chart.  The  effects  of  the  subdomlnaot  roots  are  promi¬ 
nent  in  this  region)  and,  even  for  a  linear  system,  accurate  solutions 
could  not  be  eocpected.  The  erratic  behavior  of  q^^  and  I3,  in  particular 

the  swinge  of  I3  into  positive  regions,  appears  to  lead  to  a  change  in 

the  altitude  of  the  airplane  which  results  in  a  dive  and  crash  since  no 
corrective  control  action  is  taken. 


As  a  whole  the  flight*  with  all  of  tbs  quadrature  intervals  are 
satisfactory  except  those  with  the  lOO-mseo.  Interval. 
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9*  Immolaan 

For  the  Imelman-Tam  maneuver  p  and  ifere  selected  as  the 

variables  of  primary  Intsreati  therefore^  evaluation  of  the  data  obtidxted 
ifSs  directed  toward  p  and  q^*  ^be  results  ere  displayed  in  Figures 

Ihe  results  for  p  kith  Quadrature  intervals  of  10,  $0,  and 
100  msec,  are  displayed  Ih  the  figures*  The  20-  and  70-«iseo.  results  Sxe 
not  displayed,  but  they  are  in  close  agreement  with  those  with  $0  msec* 
The  results  with  100  msec*  deviate  f  rom  the  others  during  the  period  3U 
to  1|2  seconds*  Note  that  p  goes  down  to  -0*3  at  t  ■  3I4  seconds  with  the 
O^y  The  deviation  is  less  with  the  f omuls,  again  showing  Its 

tendency  to  give  si-*-’  .I;-  better  results  for  a  violent  maneuver* 

The  results  for  q^  with  quadrature  intervals  of  10,  -and 

70  nsec*  are  displayed  In  the  figures*  The  msults  with  20  nsec*  are 
near  those  with  $0  nsec.)  and  the  results  nlth  all  four  Intervals  are  In 
close  agmsment  —  the  agreement  heing  a  little  better  with  the  0.  to°31 

than  with  the  0^^*  The  results  with  100  nsec,  are  too  large  to  show* 

They  osclUate  wildly  and  reach  both  of  the  Unit  stoppers  with  O^^and 

the  positive  limit  stopper  with  0^  toSi  during  the  early  portion  of  the 

Bianeuver*  The  subsequent  decrease  In  amplitude  Indloatec  stability^ 
howerrer* 


IV,  CONCLUSIONS  AND  RECOMMENDATIONS 


In  connection  with  the  evaluation  of  tha  accuracy  of  the 
mathematical  procedures  used  In  digital  airplane  flight  simulation^ 
selected  maneuvers  were  run  on  ODOPTT  by  means  of  an  autopilot  program. 
During  the  execution  of  a  maneuver;  digltcd  prlnt-oute  and  analog  trac¬ 
ings  were  obtained  for  the  variables  of  primary  interest.  The  digital 
samplings  were  plotted  with  the  same  time  icald  so  that  direct  compari¬ 
son  of  tlie  simulation  with  different  quadrature  intervals;  or  quadrature 
formulas;  could  be  made. 

The  results  obtained  with  some  of  the  maneuvers  indicate  that 
the  accuracy  of  the  solutions  obtained  was  being  limited  by  underflow. 

It  appears  that  some  of  the  problem  was  caused  by  overly  conservative 
scaling  in  the  F-1CX}A  program.  In  particular  it  was  noticed  that  uzider** 
flow  made  it  difficult  to  perform  a  phugold  oscillation  and  caused  the 
steady-flight  maneuver  to  have  essentially  no  variation  of  "'fll^t  condi¬ 
tions.  This  made  it  difficult  to  obtain  information  about  the  Integra- 
ticn  formulas  themaelvea  with  these  maneuvers;  particularly  at  the  smaller 
quadrature  intervale. 

As  a  whole  the  results  of  the  testa  are  very  good.  The  close 
agreement  between  the  sets  of  results  for  each  maneuver  indicates 
accuracy  in  the  solutions  by  both  0^^  and  mod  Gurk  quadrature 

formulas  at  quadrature  intervals  of  10  msec . ;  20  msec . ,  5^  msec . ;  and 
70  msec.  The  stability  and  good  accuracy  of  these  results  bad  been 
predicted  by  the  stability-chart  criterion.  The  poorer  results  at  a 
quadrature  interval  of  100  msec,  had  also  been  predicted  by  the 
stability-chart  criterion.  However,  the  accuracy  or.d  stability  obtained, 
particularly  at  h  »  7^  msec.,  was  somewhat  better  than  was  expected. 

The  0^^  mod  Gurk  quadrature  formula  gave  stable  results  for  all 
maneuvers  at  alx^quadrature  intervals,  with  the  slmolatlons  for  all  inter¬ 
vals  except  100  msec,  being  acceptable  and  reasonably  accurate.  The 
resvilts  with  h  »  100  msec,  gave  leorge,  but  stable,  oscillatlonB  in  q^ 

and  1^  for  some  maneuvers.  These  oscillations  usually  resulted  in  a  poor 
simulation  or  a  crash.  Some  of  the  wildness  in  the  solutions  with  h  %  100 
msec,  can  be  attributed  to  th^s  effects  of  the  subdonlnant  roots  of  whose 
presence  are  indicated  by  the  pro:d.mlty  of  ph  to  the  branch  contour  with 
this  quatoature  interval.  T>ie  increase  in  frequency  can  also  be  attribu*ted 
in  part  to  the  increase  in  |ih;  however  it  was  noticed  that  the  increase 
in  frequency  was  greater  than  was  indicated  by  the  stability  chart. 

The  mod  Gurk  quadrature  form.ula  also  gave  good  results 

for  all  quadrat^e'^intervals  except,  in  some  cases,  h  «  100  msec.  The 
^30^31  ^33  close  in  most  cases  giving  further 

indications  of  accuracy.  However,  the  formula  seems  to  give  some¬ 

what  better  results  with  h  «■  100  msec.^  partioularly  during  the  violent 
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portions  of  a  nanouver  and  Ml>f  the  high  fkvinenoy  pitch  dsclllations* 

This  indicates  that  giv^B  senewhat  better  resultA  for  violexit 

maneuvers* 

It  had  been  eocpeoted  that  the  numerical  Oolh^ipns  for  atrong^y 
nonlinear  systems  iioald  differ  from  the  solutions  pMUjbted  by  the  use 
of  stablllty-bhart  criteria  since  stability  charts  wer^  developed  hslag 
linear  theory*  (Refer  to  Section  II'  ••  2.  partloulahly  the  disdusslon 
relating  to  solution  by  Chj  mod  Ourk.)  Thus  the  good,  and  predictable, 
results  obtained  in  the  majority  of  maneuvers  indicates  that  both  of  the 
foonulas  used  are  not  affected  to  any  great  degree  by  nonllnearltles  of 
the  type  encountered  in  the  F-IOOA  equations.  However,  as  was  pointed 
out  in  the  discussion  of  the  Qsoillatione-in-Pltch  maneuver,  the  fre¬ 
quency  of  the  results  with  the  formula  at  h  •  100  nsec,  was  pre¬ 

dictable  by  stability-chart  oriterle*  whereas  the  frequency  of  the  results 
with  the  O33  formula  at  h  -  100  msec,  was  larger  than  the  predicted  value* 
The  results  with  the  0^^  formula  at.  h  «  100  msec*  were  also  much  larger 

in  amplitude  than  had  been  eocpeoted.  Thle  Indloatee  that  0^^  mod  Gurk 

fCrtmila  ia  affected  by  the  noinllnearltlas  in  the  F-lOQA  simulation 
eqnsiloBB  if  a  quadrature  interval  is  used  idilch  brings  the  **dimension- 
less  frequencies"  near  ths  branch  contour.  The  formula  does  not 

seem  to  be  a^grtatly  affected  by  the  nonllnearitiae,  for  the  results  mdre 
still  predictable. 

The  results  obtalne4  using  the  0^^  and  the  ^30^31 

molas  seen  to  show  that  when  formulas  of  this  type  are  applied  to  slightly 
nonlinear  systeme,  such  as  th»  F-IOOA  flight  equations,  the  stability- 
chart  criteria  hold*  This  includes  the  linitl^  condition  that  where  the 
effects  of  the  arabdomlnant  mote  are  expected  to  beodse  elgnlflcant  the 
nonlinearities  will  also  begin  to  affect  the  solution* 

The  experiments  performed,  particularly  the  Fhugold-Oseillatioa 
maneuver,  show  that  the  0^  mod  Qurk  quadrature  formula  gives  good  results 

for  systems  whobe  dimensic^ess  frequencies  lie  near  the  origin'  of  the 
stability  chart*  This  wSA  expected  frem  stablllty-ohart  criteria,  but  It 
is  of  particular  interest  because  of  the  dffleulty  oftar  encountered 
whan  using  analog  devices  for  such  systems.  The  experdnetTbs  performed, 
particularly  the  OBcillations-in-f(j.teh  maneuvers;  and  the  OsclUatiens-in- 
Isw  maneuvers,  also  demonetrated  conclusively  that  both  quadrature  formulas 
used  give  good  results  for  systems  with  frequencies  as  high  as  one  cycle 
per  second  (i.e*,  periods  as  nail  as  1  second)  with  quadrature  intervals 
as  largs  as  50  msec.  It  was  also  shown  conclusively  that  for  eystems  with 
periods  of  two  seconds,  such  as  were  preaant  in  the  OsoHlations-ln-lkv 
auneuver,  good  simulstlcns  will  he  obtained  with  quedrature  intervals  air 
large  as  100  Jiseo.  This  agrees  with  the  bshacvlor  pradleted  by  ths  sta¬ 
bility  charts. 
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Th«  prwiow  of  underflow  in  eone  of  the  menetoTere  (iH  partly 
eular  the  elwr  damonatretlon  of  Ita  pivaenoe  and  effects  In  the 
phugoitl  naneu'rer)  brings  up  the  questldn  of  iufHoient  preoisldn  in 
TnXFTT*  It  should  be  noted  thnt  tte  esdstlng  P-lOQi  sinuletlon  progran 
Is  conservatively  sealed*  The  scaling  in  the  program  is  such  that  the 
troncetiott  feature  built  into  the  nacfaine  as  protection  against  over- 
flow  is  not  used*  Until  tests  are  pezfoxmed  idth  a  less  conservatively 
sealed  progran*  it  Is  difficult  to  say  whsther  or  not  the  appearance  cf 
underfloif  indicates  a  need  for  more  precision  in  future  siaulators  of 
this  type*  The  prohlesis  Involved  in  sealing  the  progran  to  alloir  both 
violent  nsneuvers  and  aeourate  phugoids  would  indioate^  howaver*  that 
the  use  of  floatlng*point  arithsMtie  in  future  siaulators  would  be  sn 
advantage* 
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VI.  Appmbiasa 
E.  Qtiadratvire  Pottnulae 


The  B°id  Gudc  quadratulrt  fonUl*  li 

t  3  II 


^n 

*^30 

It 

3 

^31 

II 

V2^ 

II  ^  1 

and  X  ,  are  the  three  paat  ordlnatea;  k  is 
n—3  n 

ttae  preaen-^  tliae  derivative  computed  ty  ualng  tbe  value  obtained  vith 
the  O^Q  (.\)i  is  ilie  quadrature  Interval;  and 

«  I.8O65805 

-  -1.1093273 

-  0.3027468 
a  ,  »  1.1462064 
a^g  «  -0.2010870 

-  0.0548787 
b^Q  -  0.9086705 


For  uae  in  the  computer  program  the  O-aC-.^  ®od 
rature  formula  la  vrltten  ^  ^ 


where 


II 


2^*0  Vl*^  W2^  ®oV3^ 

II  It  f1 

2(a„x  .  +  b  O  + 

C  D"*!  C  11*2  c  c 


• 


a  « 


'b.  - 


« 


1-0. 

2  *o2 

3lr  A  ^ 

2  *o3 

Iv 

I  ^2 

I  “03 

I 


Gurl:  quad* 
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Thus,  is  the  only  coefficient  which  has  to  be  changed  when  the 

quadrature  interval  is  changed.  The  values  of  d^  for  the  various 
quadrature  intervals  are  given  in  Table  1. 


Quad.  Interval 

dc 

10  msec. 

20  msec. 

50  msec. 

70  msec. 

100  msec. 

0022470 

0045160 

0135030 

0202210 

0272060 

Table  1  Values  of  d^  for  the  various  quad¬ 

rature  intervals  (values  are  in 
octal  and  the  units  point  is  at  the 
extreme  left) . 
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H,  ReBiilts 


This  appendix  contains  the  tracings  for  the  results  discussed 
In  Section  III,  The  values  for  the  quantities  are  plotted  in  octal  scaled 
as  obtained  from  the  con9>uter.  The  notation  Bn  (n  •  0,  1,  2,  20) 

indicates  that  the  number  shown  must  be  multiplied  by  20  to  obtain  the 
true  value. 
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